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ABSTRACT 


Space Shuttle Orbiter landings indicate both long term directional instability and the 
potential for pilot induced oscillations during landing and rollout before nosewhee! touchdown. 
The Orbiter’s Flight Control System requires improvements to increase directional contro! in the 
two point stance (after main gear touchdown with the nose in the air). A number of modifications 
are proposed to improve directional control. This thesis describes the control deficiency, potential 
improvements to the Flight Control System (FCS), and evaluates a number of these 
improvements. The evaluation was performed by modeling the Orbiter’s postlanding 
lateral/directional control laws using a commercially available engineenng software package 
known as MATLAB 5.0. Directional control of the Orbiter was evaluated with and without the 
proposed modification to obtain a comparison of control response. 

Initial evaluation of future Orbiter FCS modifications could be performed using 
commercially available engineering software packages such as MATLAB; rather than costly full- 


up Orbiter simulators. A low cost initial evaluation of changes may save NASA resources. 
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I. INTRODUCTION 


The Space Shuttle, operated by the National Aeronautics and Space Administration 
(NASA), is the only manned space flight vehicle currently operated by the United States. The 
Space Shuttle consists of the Orbiter, External Tank (ET), and two Solid Rocket Boosters 
(SRB’s). The Orbiter contains the astronauts and payloads. It is the only piece of the Space 
Shuttle that completes the entire mission. The Orbiter returns from space to land similar to a 
conventional gliding airplane. The continued safe operation of the Space Shuttle is the number 
one concer of the Space Shuttle Program (SSP). A recent safety concem its the directional 
control of the Orbiter during one segment of the landing/rollout phase [Ref. 1]. This perceived 


directional control deficiency is described below. 
A. BACKGROUND 


During the landing and rollout phase of flight, the gains associated with the piloting task 
are greatly increased. The margin for error is vastly decreased as pitch, roll, and yaw requirements 
become more stringent through the flare-to-touchdown and subsequent rollout. During this 
critical phase, pilot control commands will be more frequent and often of greater magnitude to 
make timely corrections to flight path and counter external perturbations such as crosswinds, wind 
shear, or gusty conditions. These increased demands on the control system can occasionally result 
in a condition known as Pilot Induced Oscillation (PIO). 

PIO occurs when a pilot attempts a particular control task which is beyond the capabilities 
of the aircraft control system. PIO can only occur with the pilot actively in the vehicle control 
loop and is therefore not a stability concern. A PIO is typically manifested by an uncommanded 
oscillatory motion that can only be stopped when the pilot reduces control inputs or “backs out of 
the loop.’ The severity of PIO can range from an annoyance to complete loss of control. The 
result of a PIO situation can often be unpredictable. What might be an annoyance dunng landing 
one day, might be disastrous on the next. PIO can result from many different root causes 
including an excess amount of friction in a mechanical control system, actuator rate limiting, or 
transport delays in today’s most advanced flight control systems. In all cases, a large phase lag 
from the control input to aircraft response, lies at the heart of the problem. 

The Orbiter incorporates a ‘fly-by-wire’ control system to command the flight control 


surfaces, depicted in Figure 1. 


Inboard Elevon 


Body Flap 


Rudder/Speedbrake 


Outboard Elevon 
Figure 1. Orbiter Flight Control Surfaces 


The flight controls, control stick, rudder pedals, and speed brake switch are not 
mechanically connected to their respective control surfaces. Orbiter control 1s afforded through a 
myriad of feedback control loops, filters, gain multipliers, integrators, differentiators, etc. The 
goal of these components, in conjunction with a feel system, is to supply the pilot with an Orbiter 
which responds acceptably to all control inputs. Fly-by-wire control systems, if not carefully 
designed, can be susceptible to PIO. 

The directional control problem of interest is most evident after touchdown while the 
Orbiter is in the two point stance with the mainwheels on the ground and the nosewheel in the air. 
This flight segment, from initial touchdown until derotation is complete, 1s approximately 20 
seconds in duration. With dispersed touchdown speeds between 185 and 225 knots, a small 
deviation in track can result in a runway departure very quickly. Recent Orbiter landings have 
shown that there is a directional PIO problem in the two point stance[Ref. 1]. The problem seems 
to abate after nose gear touchdown due to the additional directional control afforded by the Nose 
Wheel Steering (NWS) system. The PIO is associated with the task of capturing and tracking the 
runway centerline and may be exaggerated in conditions of greater crosswind. The possibility of 
the Orbiter becoming directionally uncontrollable and departing the runway is unacceptable. This 


problem raises some critical questions including: 


1. Are there operational impacts? Must the crosswind limit be reduced to protect against 
PIO? Should Shuttle landings require wider ninways and/or lakebed landings? 

2. Are there hardware modifications such as changes to the rudder actuators that would 
improve directional control? 

3. Are there potential software modifications to the control laws which will reduce or 
eliminate this PIO condition? 

4. If PIO is evident in a nominal landing, what would happen in the case of a blown 


main gear tire? 
iP Operational Impacts 


Reducing the crosswind limit may reduce the potential for direction PIO; however, it will 
not eliminate the problem. In fact, the landing which experienced the most significant directional 
PIO had a 2 knot crosswind [Ref. 2]. Other flight data also shows that even with minimal 
crosswind, directional PIO can occur. 

In addition, there is appreciable pressure to increase the crosswind limits for the following 
reason. The commitment of the United States to building the International Space Station has 
extended the mission of the SSP from pure science to upholding national political obligation. The 
proposed International Space Station will be in a high inclination orbit which will reduce launch 
opportunities and decrease the length of launch windows. Landing weather conditions are a 
concem on launch day since the potential for an emergency landing such as a Return To Launch 
Site (RTLS) or Trans-Atlantic Landing (TAL) can occur on launch day. Meeting the Space 
Station construction schedule would be much easier if launches and landings could be performed 
in a greater variety of wind conditions. 

Requiring lakebed landings for all Shuttle flights 1s unacceptable. Orbiter landing mass 
properties are greatly restricted for lakebed landings to maintain nose gear loads within limits 
[Ref. 3]. Most Orbiter landing mass properties require concrete runways. Further, landing on a 
lakebed requires landing at either Edwards Air Force Base (EAFB) in Califorma or Northrup, 
New Mexico which 1s impossible for most ascent aborts modes. 

The Shuttle Landing Facility (SLF) at Kennedy Space Center (KSC) and the concrete 
runway used for Shuttle landings at EAFB are 300 ft. wide. The SLF has 50 ft. of usable shoulder 
on each side. In comparison, some contingency and abort runways are less than 150 ft. wide. 


Obviously, it is impossible for the SSP to widen all runways available for landing. Restricting the 


Orbiter from landing on runways considered too narrow, greatly reduces abort coverage (no where 


to landing with an SSME failed) to an unacceptable level. 
pep Hardware Modifications 


Hardware modifications, such as redesign of the rudder actuator to increase deflection 
rate, may reduce susceptibility to PIO. However, hardware modifications are expensive, require 


recertification, and demand vehicle down time to implement the changes. 
3. Software Modifications 


The obvious long term solution to the PIO problem lies in redesigning the existing flight 
control laws to afford better directional control in a wider variety of wind conditions. 

A redesign of this magnitude would traditionally require the control engineer to estimate 
the new control laws followed by experimental Verification and Validation (V&V) of the new 
control laws in the Vertical Motion Simulator (VMS) located at the NASA Ames Research Center 
(ARC) located in Sunnyvale, California. This procedure typically entails qualified astronauts 
traveling from the Johnson Space Center (JSC) in Houston, Texas to the ARC VMS to fly 
numerous simulated approaches in varying wind conditions. This evolution can be time 
consuming and costly. 

With the advent of faster more powerful personal computers and workstations, and the 
availability of very capable software packages which can operate in both the frequency and time 
domains, the possibility of wnting and evaluating new control laws immediately has become a 
reality. This breakthrough can eliminate the costly iterative process of tweaking and reevaluating 
control laws in an expensive motion base simulator. The final V&V is still necessary, but the total 
time and cost are drastically reduced. 

The use of simulation to evaluate control laws rely heavily on the numerous models used 
to develop the simulation. A normal flight test program would not rely solely on simulation 
results for solutions to PIO problems. The Shuttle Program must rely on these simulations results 
due to the nature, cost, and risk of flying Shuttle flights only to evaluate control law changes. 
Flight data is gathered from each flight; however, flight conditions do not test the extreme ends of 


the envelope. There is no possibility of performing a complete flight test program. 


PURPOSE 


This thesis has two objectives: 

(1) Develop changes to the control laws to reduce the directional PIO tendency present 
between mainwheel touchdown and completion of derotation while the Orbiter is in 
the two point attitude. 

(2) Evaluate these changes using a Commercial Off-The-Shelf (COTS) PC-based 
software package. 





Hl. ORBITER DIRECTIONAL STABILITY AND CONTROL 


The nominal Orbiter rollout profile, directional stability, and directional control in the two 


point stance are described in this chapter. 
A. NOMINAL ROLLOUT PROFILE 


The nominal touchdown and rollout profile, with approximate speeds and time from 


touchdown, is described in Table I[{Ref. 2, 3 and 4]. 


Table 1. Nominal Touchdown/Rollout Profile 


[Fochaoun” moss [0 
[Binks fncons™ | a0 | 


Commence auto derotation™ 11.0 
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Commence manual 175 16.0 
derotation™ 
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Drag chute fully disreefed | = 175, 16.0 
Nosewheel touchdown 145 to 155 19.0 


First brake application 140 (KGS) 23.0 


Drag chute Jettison‘”’ 0 +20 (KGS) 34.0 


Notes: 1. Targeted touchdown speed of 205 KEAS for Orbiter gross weights 2 
220K Ibs. (195 KEAS for Orbiter weights < 220K Ibs.). 
2. Built In Test (BIT) requires 1.9 sec. after Weight On Wheels Lock-On 
(WOWLON). 





3. Chute deploy after main gear touchdown. 

4. Nominal derotation is automatic with beep tnm activation targeting 1.5 
deg/sec. 

5. Only performed if beep tnm fails. 

6. Brake initiation at 140 KGS or 5,000 ft. runway remaining, whichever occurs 


first. 


7. If drag chute is not jettisoned by 40 KGS, the chute will remain on to 


minimize damage to main engine bells. 
B. DIRECTIONAL STABILITY IN THE TWO POINT ATTITUDE 


Directional stability in the two point stance is different than the directional stability of the 
Orbiter in flight. In flight, the Orbiter would typically rotate about the Center Of Gravity (COG). 
In the two point stance, assuming no lateral skidding of the mainwheel tires, the Orbiter 1s 
constrained to rotate directionally about a point on a line that runs laterally through the areas of 
contact between the mainwheels and the runway. This point is referred to as the Center Of 
Rotation (COR). Directional stability is governed by four main factors in the two point attitude: 

(1) The relative position of the COG to the COR. 

(2) The increased friction on the downwind mainwheel caused by the rolling moment 
generated as a result of the dihedral effect or roll due to sideslip, Lg, generated 
primarily by the highly swept wing. 

(3) The aerodynamic ‘weathercock’ effect or yaw due to sideslip, Ng, generated primarily 
by the large vertical tail. 

(4) The relative position of the drag chute to the COR. 

The first factor will attempt to keep the Orbiter traveling in a straight line with respect to 
the runway in a zero crosswind condition but can act to direct the Orbiter away from a crosswind 
as a result of the second factor. The third factor will tend to fair the Orbiter into a crosswind in 
flight. After touchdown, with the COR shifting aft, this ‘weathercock’ effect becomes a much 


weaker contributor. The forth factor will attempt to fair the Orbiter into a crosswind. 
1. Center of Gravity (COG) Effects 


The nominal COG envelope is 32 in. long and 3 in. wide, centered laterally and 
positioned approximately 8 4 ft. forward of the mainwheels [Ref. 5]. The COG will always be 
well ahead of the mainwheels in the two point stance. The COG geometry with respect to the 


mainwheels is depicted in Figure 2. 
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Figure 2. Orbiter Center of Gravity Location 


With this geometry, the COG acts as a pendulum in the honzontal plane. The 
deceleration of the Orbiter caused by the friction in the mainwheels and the mainwheel brakes acts 
on the COG to produce a force in the Orbiter’s direction of motion. This force is analogous to the 
gravitational force on a normal, vertically onented, pendulum. The result is a restonng yaw 
moment which is always attempting to maintain the Orbiter on a straight track in the direction 
from the COR through the COG. In a zero crosswind condition with no mainwheel braking, the 
COR will most likely be positioned directly between the mainwheels. In this event, the Orbiter 
will have positive directional stability in the direction of motion. In a crosswind case, the high 


dihedral of the Orbiter can become a directionally destabilizing factor. 
Zs Dihedral Effect 


The Orbiter has a large dihedral effect due primarily to the highly swept wing. Dihedral is 
the tendency for an aircraft to roll away from sideslip. In a landing scenario, the sideslip angle can 


be calculated as follows: 


crosswind component 


=— ta =: nn 
i cties: Orbiter Indicated Airspeed 


The rolling moment generated in a highly swept wing aircraft in response to sideslip can 
be explained by two factors. First, the upwind wing senses flow closer to perpendicular to the 
leading edge which decreases flow in the spanwise direction; thereby increasing lift. The 
downwind wing has the opposite effect of increased spanwise flow; therefore decreased lift. The 
second, compounding effect, is a blanking effect on the downwind wing caused by the large 
fuselage structure of the Orbiter. This reduction in airflow also decreases lift. 

Another contributor to positive dihedral is the high vertical tail. The center of pressure 
associated with the airflow across the vertical tail is well above the COR which would cause a 
further roll away from crosswind. The net result is a rolling moment away from crosswind as 


depicted in Figure 3. 
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Figure 3. Orbiter Dihedral 


This rolling moment will cause an increased down force on the downwind strut and a 
lighter force on the upwind strut. This differential load will cause increased friction on the 


downwind wheel and a subsequent shift in the directional COR to the downwind side of Orbiter 


centerline. This shift in the COR, in conjunction with the COG effect discussed in subsection 1, 


will act to cause the Orbiter to deviate away from crosswind as depicted in Figure 4. 
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Figure 4. Orbiter Yaw Due to Dihedral 


3. Weathercock Effect 


The weathercock effect is typically a tradeoff between the destabilizing effect of the 
fuselage and the stabilizing effect of the vertical tail. The aft location of the COG of the Orbiter 
in flight is somewhat destabilizing since the stabilizing moment arm associated with the fin is 
shorter and the destabilizing moment arm of the fuselage is longer. When the Orbiter attains the 
two point stance, the COR shifts even further aft. This change creates a condition of weak 


instability. The Orbiter will tend to point away from a crosswind as depicted in Figure 5. 
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Figure 5. Orbiter Yaw Due to Weathercock 


4. Drag Chute Effect 


The drag chute was incorporated on to the Orbiter to decrease landing rollout distances for 
certain contingency runways. Along with the decreased rollout, there is a small amount of 
positive weathercock effect. It 1s important to note that the drag chute is nominally fully disreefed 
three seconds pnior to nosewheel touchdown as described in section II.A [Ref. 2 and 3]. 
Therefore, the drag chute will only be a brief factor to steady-state directional stability in the two 
point stance. Yet the three second impulse may be enough to perturb the steady case. The drag 
chute will blow downwind from a crosswind. Since the drag chute is attached aft of the center of 


rotation, the resulting yaw moment on the Orbiter is into the crosswind as depicted tn Figure 6. 
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Figure 6. Orbiter Yaw Due to Drag Chute 


Si Frequency and Damping 


A summary of the directional stability parameters is shown in Table 2 along with their 


relative strengths. 


Table 2. Summary of Directional Stability Parameters 


Dihedral Away from crosswind 


Drei | Aarne [Weak [encom 


All of these factors combine to yield fairly weak directional stability in the two point 











stance. This weak stability, when modeled with a simple second order system, would yield a low 
open loop frequency response to a disturbance. While this data point has never been collected in 
flight, a rough estimate of what one might call ‘the directional natural frequency’, ®p, can be made 
by observing video of actual closed loop cases of PIO. The videoThe open loop natural 


frequency, Wp, is estimated to be on the order of | rad/sec. The only yaw damping, N,, 1s most 
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likely due to the large vertical tail and possibly any differential friction in the mainwheels. Both 
of these factors would be small which yields very low directional damping. The complete 
dynamic stability of the Orbiter in the two point stance can be modeled by the following second 


order equation where Np is total damping and Ns is total stability. 


8 + Nast No 
where Wp = yf N73 


ae Nip 
and the damping ratio, ¢, = 


Zao, 





The natural frequency and damping ratio will be estimated from the simulation results. A 
low frequency and low damping ratio will indicate the tendency of the FCS to be susceptible to 
PIO. 


C. DIRECTIONAL CONTROL IN THE TWO POINT ATTITUDE 


Directional control in the two point stance can be affected with two controllers; 
aerodynamically with the rudder, and mechanically with differential braking. Brakes are not 
employed until derotation is complete to avoid excessive nosewheel slapdown rates and to reduce 
the risk of a blown main tire. Braking nominally begins at 140 knots [Ref. 3]. Differential 
braking would only be employed in extreme situations, such as a flat main gear tire. Therefore, 


directional control in the two point stance would be afforded only by the rudder. 
i Rudder 


The rudder is actually a dual purpose control surface acting as both a directional controller 
and a speedbrake as depicted in Figure 7. 
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Figure 7. Orbiter Rudder/Speedbrake 


The surface is mechanically limited to a maximum deflection rate of 16 deg/sec normally 
[Ref. 6]. Postlanding the deflection rate is further limited tol 3deg/sec since the speedbrake is full 
open. The software limits the deflection rate to 14 deg/sec so that sufficient hydraulic fluid is 
available to other control surfaces. In the event of the lost of one or two hydraulic system, the rate 
is limited to 12 deg/sec and 7 deg/sec, respectively [Ref. 4]. The surface areas is approximately 
97 ft. and can generate up to 2.25 lateral g’s. 


ie Directional FCS Susceptibility to PIO 


The FCS has many limitations which result in the high potential for directional PIO in the 
two point stance. The slow surface rate of the rudder, slow response time, high control power of 
the rudder due to the large surface area, and rudder placement with respect to the gear are 
conducive to PIO. If a large pedal input is sustained until the yaw rate is well underway, loss of 
control is probable. The slow surface rate does not allow an opposite input to be effected quicky 
enough to counteract the yaw rate. 

In addition, ARC VMS studies have indicated rudder rate saturation is encountered in 


high crosswind landings or in performing postlanding lateral maneuvers stress tasks [Ref. 7]. 





Il. RUDDER CONTROL 


This chapter describes the rudder control path including the software command path and 
rudder actuator. Detailed drawings and descnption of the Orbiter FCS control laws are found in 


the Software Requirements Document [Ref. 8] and Lockheed drawings [Ref. 9]. 
A. RUDDER COMMAND PATH 


The control logic for the rudder surface may be broken down into two basic steps. The 
first converts rudder pedal deflection angle into lateral acceleration error command. The Orbiter’s 
flight control system uses lateral acceleration feedback to null errors to the commanded lateral 
acceleration and provide directional stability augmentation. The second step converts the lateral 
acceleration error command into a yaw rate which is filtered of Orbiter structural bending modes 
and finally converted into a command to the control surface for degrees of deflection. This 
command is sent to Prionty Rate Limiting (PRL), which limits the rudder deflection rate 
depending on the number of hydraulic systems available. These two steps will be henceforth 


referred to as the rudder command feedback and rudder deflection command portions. 
I. Rudder Command Feedback 


A simplified block diagram of rudder pedal position conversion to lateral acceleration 


error command is shown in Figure 8. 
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Figure 8. Rudder Command Feedback 


The various components of the rudder command feedback path are discussed in the 


subsections below. 


a. Rudder Pedals 


There is a Rudder Pedal Transducer Assembly (RPTA) located at both the 
Commander and Pilot stations. Each RPTA consists of two pedals and a transducer which 
converts degrees of rudder pedal deflection into electronic signal. The RPTA uses a time constant 
of 100 ms in the mechanical to electronic transfer. This signal is processed at 6.25 Hz and passed 


to the rudder command shaper with the name tag DRMAN. 
b. Rudder Command Shaper 


The rudder command shaper converts DRMAN into a shaped and scaled rudder 


pedal command named DRMS as shown in Figure 9. 
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Figure 9. Rudder Command Shaper 


The shaping schedule is parabolic and conforms to the following function: 
DRMS = (0.042)DRMAN? + (0.131) DRMAN . A maximum deflection limit of 22.5 degrees 


is imposed on the output signal. 
cS Rudder Pedal Gain Scheduler 


The rudder pedal gain scheduler produces a gain signal, which when multiplied 
with the shaped and scaled rudder pedal command, DRMS, produces a lateral acceleration 
command. This gain is scheduled as a function of True Airspeed (TAS) and is given the name 
GNYDRM as depicted in Figure 10. This gain also converts the units of the signal from degrees 
of pedal deflection to lateral g’s. The rudder pedal gain scheduler samples at 1.04 Hz. 





Figure 10. Rudder Pedal Gain Scheduler 


a. Rudder Trim Command 


The rudder trim command signal is generated by converting the command signal 
from the rudder trim panel into a commanded lateral acceleration bias. This is performed with a 
rectangular integrator which ramps the output signal and also limits the maximum bias to 0.1 g’s. 
This signal is sampled at 6.25 Hz. The output signal is the lateral acceleration trim signal, DRTI, 
which is summed with GNYDRM from the rudder pedals to produce a total lateral acceleration 
command, DRTMS, as depicted in Figure 11. 
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Figure 11. Rudder Trim Commands 


é. Lateral Accelerometers 


There are four Accelerometer Assemblies (AA) located forward of the crew 
compartment. These AA’s directly measure lateral acceleration and simultaneously filter the 
signal as shown below: 


Measured Acceleration 
—> Lateral Acceleration Signal( NYP) 





(rata same 
s* +2(.5)(15)s+15° /\s? + 2(.8)(146)s + 1467 


This signal is sampled at 25 Hz and passed to the lateral acceleration feedback 


compensator. 
fs Lateral Acceleration Feedback Compensator 


The lateral acceleration feedback compensator compares the actual lateral 
acceleration, NYP, as passed from the AA’s, to the commanded lateral acceleration, DRTMS, to 
generate an error signal, or additional lateral acceleration command. A low pass filter is 
incorporated in this procedure to smooth the output. The resulting signal is sampled at 25 Hz. 
The process is depicted in Figure 12. 
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Figure 12. Lateral Acceleration Feedback Compensator 


2: Rudder Deflection Command 


A simplified block diagram of lateral acceleration conversion to rudder deflection 


command is shown in Figure 13. 
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Figure 13. Rudder Deflection Command 


The various components of the mdder deflection command path are discussed in the 


subsections below. 
a. Lateral Acceleration Gain 


The lateral acceleration gain scheduler produces a gain signal, GRAY, which 
when multiplied with the lateral acceleration error command, DRTMS, produces a yaw rate 
command, DAY. This gain is scheduled as a function of Mach as depicted in Figure 14. This 
gain also converts the units of the signal from g’s to deg/sec. The sample rate for this scheduler 1s 


1.04 Hz. 





Figure 14. Lateral Acceleration Gain Schedule 


b. Cosine of the Angle of Attack (a) 


A dimensionless signal, cos(a), is mixed with the yaw rate command signal and 
summed with the yaw rate error to produce a yaw rate command about an inertial axis. This 


signal will always reduce the yaw rate command signal in producing a yaw rate with respect to the 


runway. The cos(a) is sampled at 1.04 Hz. The output signal from this mix is referred to as the 


unfiltered yaw rate command, DRRCUF. 
Cc. Rate Gyro Assemblies 


There are four Rate Gyro Assemblies (RGA’s) located on the aft bulkhead below 
the floor of the payload bay. Each of the RGA’s sense rotation about the three vehicle axes and 


simultaneously filter the signal as shown below: 
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These rates are the primary feedback to the flight contre] system to maintain 


control. 
d. Yaw Rate Structural Filter 
A second-order structural bending filter 1s applied to the yaw rate feedback. 
é. Yaw Rate Error Determination 


The filtered body yaw rate, DRPRM, is summed with the product of the sideslip 
angle times the cosine of the angle of attack, BCOSALF, and yaw rate command, DAXFDC, to 
determine the yaw rate error, DRERR, as shown in Figure 15. The sideslip angle is not used in 
the approach and landing flight phase. Therefore, BCOSALF is zero in this analysis. The 
resulting signal is sampled at 25 Hz. 
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Figure 15. Yaw Rate Error Determination 
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ie Structural Bending Filter 


A structural bending filter is included in series to remove any known Orbiter body 
bending modes. The resulting signal is the yaw command, DRRC, which is processed at 25 Hz. 
The filter is made of two consecutive second order notch filters. This fourth order system will 


inherently result in transport delays which would only contnbute to the time lag problem. 
g. Yaw Rate Gain 


The yaw rate gain scheduler produces a gain signal, GDRC. This gain and a 
constant, GERYFWDS, 1s multiplied with the yaw command signal, DRRC, to produce an initial 
rudder deflection command, DRCPF. This gain is scheduled as a function of Mach and dynamic 
pressure as depicted in Figure 16. This gain also converts the units of the signal from deg/sec to 


deg. The sample rate for this scheduler is 1.04 Hz. 
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Figure 16. Yaw Rate Gain Scheduler 


h. Rudder Deflection Limiter 


The rudder deflection limiter provides a maximum rudder deflection limit in the 
positive and negative directions. The deflection limit, DRC LIM, is 24.1° and 27.1° in flight and 
post Weight On Wheels Lock-On (WOWLON), respectively. Within these limits, the limiter 
passes a one-to-one correlation from the filtered rudder deflection command, DRCPF, to the 
actual rudder deflection command, DRC. The WOWLON signal is sampled at 6.25 Hz, but the 


rudder deflection limiter processes at 25 Hz. The system is configured as shown in Figure 17. 
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Figure 17. Rudder Deflection Limiter 


L. Priority Rate Limiting (PRL) 


The rudder deflection command is sent to PRL. PRL limits control surface 
deflection rates to manage the amount of hydraulic fluid sent to each of the actuators. There are 
three Auxiliary Power Unit (APU)/hydraulic systems available to provide hydraulic fluid to the 
control surface actuators. Normally, with three hydraulic systems available, PRL limits the rudder 
deflection rate to 14 deg/sec [Ref. 4]. This slow deflection rate is conducive to the pilot’s rudder 
inputs being out of phase with Orbiter response. The deflection rate is further reduced depending 
upon the number of hydraulic systems failed. With one and two systems failed, the rudder 


deflection rate is limited to 12 and 7 deg/sec, respectively. 
B. RUDDER ACTUATOR 


The rudder/speedbrake actuators receive hydraulic fluid from all three hydraulic systems 
through a switching valve. The switching valve produces one hydraulic output pressure source for 
the rudder/speedbrake servovalves. The power drive unit, mechanical actuator, and actuator drive 


shaft are shown in Figure 18. 
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Figure 18. Rudder/Speedbrake Actuator Location [Ref. 6] 
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IV. EQUATIONS OF MOTION AND STABILITY AND CONTROL 


DERIVATIVES 


A. ASSUMPTIONS 


The following assumptions were used in developing the Equations of Motion (EOM) for 


the Orbiter in the two point stance: 


EF 
a 


Inertial frame fixed to the flat Earth since the Orbiter is on the ground. 

Constant vehicle mass properties since there is no propellant consumption and 
nothing ejected until drag chute deploy. 

Vehicle is symmetric about the XZ plane. Therefore, I,y = lyz = Cyo = Cio = Cro = 9. 
Equations of motion are linear about a steady state flight condition. Small 
perturbation theory implies that the squares and products of perturbations may be 
neglected. For every perturbation angle a, Sin (a) = a and Cos (a) = 1. 

Atmospheric property variations are negligible. Air density, temperature, pressure 
changes are negligible over a twenty second period and less than 10,000 feet rolling 
distance on the ground. 

Steady state rotation rates, side velocity, and bank angle are all zero. The Orbiter in 
the two point stance will not perform rapid maneuvers or fly with significant steady 
state sideslip. 

The air flow is quasi-steady. In other words, the pressure distributions adjust instantly 
to rate changes. This assumption is generally valid for angle of attack rate a if the 
following condition is satisfied [Ref. 10]: a c/2U, < 0.04. For the Orbiter, the chord 
length c = 39.56 ft. and the initial velocity U,; = 353 ft./sec. The inequality is 
satisfied if a <0.714 rad/s = 40.9 deg/s. In the two point stance prior to the initiation 
of derotation, o is essential constant. After initiation of derotation, a is 1.5 deg/sec. 
The stability derivatives with respect to « may be neglected. 

Motion in the Z direction is constrained due to the runway. Therefore, the downward 
velocity W=0 and the associate forces and moments are zero. 

Vehicle is at constant a. Therefore, all pitching moment derivatives are zero, the 


pitch rate, q, is zero, and all control derivatives with respect to q are zero. 


2] 


10. The elevons are at trail (0°) at touchdown and remain there until the initiation of 
derotation [Ref. 2]. Therefore, all variations of the coefficients due to elevator and 
aileron are zero. 


11. Forward velocity, U is constant. 
B. EQUATIONS OF MOTION 


The general EOM in the body-fixed coordinate system obtained from Roskam [Ref. 10] 
are defined below. The subscript denotes the axis and the apostrophe denotes time derivative. All 
symbols are defined in the List of Symbols, Acronyms, and/or Abbreviations on page Xi. 
Although the control problem under investigation is a directional control problem, both the 
longitudinal and lateral/directional equations are included in the derivation and simulation for 
completeness. The general force, moment, and kinematic equations are: 

F, = m(U -VR+WQ)+mgSin(©) 

F, = m(V +UR-WP)-mgCos(©)Sin(®) 
F, = m(W -UQ+VP)-mgCos(®)Cos(O) 
L = IP -IeR -IePQ+(In-ly)RQ 

M = 1,,Q +(loclz)PR+1,.(P7-R’) 

N = I,R -lgP +(lyy-lxx)PQ+LcQR 

P= @-y Sin(@) 

Q = ©Cos()+y Cos(O)Sin(®) 

R = y Cos(®)Cos()-© Sin(®) 


ele 


Each variable is the sum of a steady state components, denoted with a subscript “1”, and 
perturbation components, denoted with lower case: 
U =U,tu V=V)tv W = W,t+w 
ae Q=Qitq R=Rjt1r 
Y=Pity ©=0,+0 P= @,+o 
YP=Pitw ©=0,+0 OD=0,+6 
Peta fe 0) Ryser 
L=L)41 M=M,)+m,  N=N,+n 
The steady state forces, moment, and kinematic equations are as follows: 


F,, =m(U ,-V,;R\+W,Q))+mgSin(©)) 
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Fy, = m(V |+U,Ri-W,P))-mgCos(©,)Sin(®,) 
F,, = m(W \-U,Q,+V;P,)-mgCos(®,)Cos(®,) 
Ly = IP ploR 1-hoP1 Qt e-hy)RiQy 
My = 1,,Q 1 +(Llacla)PiRithe(P1’ -Ri’) 
Ni = LR i -hoP 1 +(lyy-lxx)P1Qi +H QiRi 
P; =@,-w Sin(©) 
Q, = ©’, Cos(®,)+y ;Cos(©,)Sin(®)) 
R, = w ,;Cos(©,)Cos(®,)-© ,Sin(®)) 
The pitch and bank angles are defined such that: 
Sin(8) = 0 Sin(o) = o Cos(8) = Cos(o) = 1.0 
This small angle assumption restricts the pitch and bank angle perturbations to less than 
15 degrees. This limit is reasonable for landing and rollout and does not constitute any serious 
restriction. Perturbations are assumed to be sufficiently small such that products and cross- 
products of the perturbations are negligible with respect to the perturbation itself: 
PV = pw = pq = pr= qu=qw=qr=ru=rv=0 
op —1 70 -Wo—-Wwio—o0 —0 
The perturbation force, moment, and kinematic equations are as follows: 
f,, = m(g8Cos(@)+W)q-V)1+u -R,v+Q,w) 
f, = m(-g@Cos(®,)Cos(@;)+g20Sin(M,)Sin(©,)-W pt+U,r+R,ut+v -P;w) 
f, = m(g@Cos(O,)Sin(®,)+g9Cos(®,)Sin(@,)+V,p-U,q-Q,u+P)v+w ) 
L= Lap (QiptPigyhot(-ly+le)(Rig+Qit)-het 
Ma = 21ePipt(Leclz)Riptlyq +(oclz)P11-21,.Rir 
n= (-lotHy)(QiptPig)-Lop +(RiqtQinletLt 
p = OF Cos(@,) - w Sin(@,)+@ 
q = 6 Cos(@,)+wCos(©,)Sin(®,)-6'¥ ,Sin(®, )Sin(O,)+¥ ;Cos(®,)Cos(©,)-O )Sin(®,))@ 
r= -(© ;Cos(®,)+'¥ ;Cos(Q,)Sin(®, ))@t+y Cos(@,)Cos(©,)-8 Sin(®)-O'¥ ;Cos(®,)Sin(@)) 


1. Longitudinal Equations of Motion 


The longitudinal forces in the stability axis system are defined as: 
Fy ee? Foy ma D, = Cpidayn> Ly = Ci iqdynS M, a Ci QaynSC 


The non-dimensional stability and control derivatives in the stability axis are defined as: 
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The longitudinal dimensional stability axis derivatives are defined as: 


Axial force: 





, = G4yS(Co, +2C, ) | é ~ 44n5(Co, +2C, | fi _ ae 

a mU, sccm m rad sec’ " OU see 
y ata, ft aS Coy, fl 
47° 2mU,  radsec 7 * m rad sec” 


Normal force: 


, Fav (Cy, +26) 1 7, = danS(Cs, #2C) fn 


; mU, sec m radsec’ ” “U;, see 





_TFa8Cr, fi, _ IanO, fl 


: 2mU, rad sec es m rad sec” 


Pitching moment: 














_ danSA Cy, +2C,, | rad ae Jan9CCm, | ae M, rad 
av eau! fisec —— sec’ "  U,_ fisec 
yee) y l 
ee GaSe Gm ae game gnSCCm, 1 
1 SoU emesee Seine SC 


Recall that motion in the Z direction is constrained due to the runway. Therefore, the 
downward velocity W=0 and the associate forces and moments X,, Zy, and M, are zero and not 
included in the state space equation. 

The pitch moment of inertia changes slightly when the Orbiter touches down. The center 


of rotation shifts from the center of mass to the axis formed by the line between the points where 
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the main gear tires meet the runway. The parallel axis theorem [Ref. 11] allows a simple shift of 
the moment of inertia with the following equation: /,,,. = /,, + mh? where lyme is the moment of 
inertia about the main gear, m is the mass of the Orbiter, and h is the distance between the center 
of mass and the new rotation point. All pitch moment dimensional stability derivative equations 
will be adjusted accordingly by multiplying by 1,,/Iyme. 


The following are the perturbed longitudinal equations of motion: 


ar at Ag gcos(@, )O+ X 5 O, 


ZL SEO} (u, +Z,)q gsin(O, )@ Z;, 6, 
st 
Ui, U, Ui, U, U 


q= M,u+ M,a+M,q+ M, 6, 








a= ¢ 

In state space form, X’=AX+BU: 
u' A, Xo xe — gcos(@, ) u Ag 
wa Z, Ly “q - gsin(@, ) peli? 

we +] fee |6 

mee UU; Up Pelee 
a M, M, ™M, 0 0 é, 

Oo 0 0 0 

ee Lateral-Directional Equations of Motion 


The lateral-directional forces are defined as: 
Fyi a ©. Gav Ly a Ci1Qdynob N, = CrGdynb 


The non-dimensional stability and control derivatives are defined as: 








6 e Oo. 
cr = 3p ey az ge Cy oo aa 3 rb C, CH 75 er ~ A. y 
2G. 2Ui 
6 CO. é é 6 
C,, = BB C, = pe C, rb C Ci, =A Ci, ~ 5, ! 
ZC 2U1 
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The lateral-directional dimensional stability and control derivatives are defined as: 





Side force: 
dan3C,, ft . q an SOC,, fi - Dams0C, sf 
Pm — radsec? P " 2mU, radsec " = 2mU,_ radsec 
dan0Cy, fi danCy, ff 
On m rad sec? a m  radsec? 


Rolling moment: 














PanS0C,, | ; Pam80°C, | Jan 50°C, 1 
a eeecce P~ 21.U, sec 9" | 2iaimeesee 
JanS6C,, | danS6C,, | 
oo  MESCSSCS@#@?? eMC” 
Yawing moment: 
Pan S0C,, , dpaS! Caml z Tiere 4 
2 sce: Pp -21,,U, sec 72 | SO sce 
damS0C,, | GamS0C,, 1 
rs re 
a ie sec? or i sec? 


ee — 


The yaw moment of inertia changes slightly when the Orbiter touches down. The center 


of rotation shifts from the center of mass to a point between the main gear. It is assumed that this 


point will be directly between the main gear. The parallel axis theorem [Ref. 11] allows a simple 


shift of the moment of inertia with the following equation: /,, = /.. + mh? where Img 1S the 


moment of inertia about the main gear, m is the mass of the Orbiter, and h is the distance between 


the center of mass and the new rotation point. All yaw moment dimensional stability derivative 


equations will be adjusted accordingly by multiplying by 1/lame. 


The roll moment of inertia is assumed to remain referenced to the center of mass as the 


main gear struts are free to compress and extend which permits the Orbiter to roil about a point 


near the center of mass. I, will not be adjusted. 
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The following are the perturbed lateral-directional equations of motion: 


ie a 
cin 


cos|Q O 
| ee 
p=L,B+t L pt L,r+—7— + b5,0, + Ls 6, 


Pa 
Y + 
U, U7 


0 


i 
ete pet N Gr a cg OnetuN 0: 
g'= p+ tan(O, )r 
For small ©’s and small I, inertial coupling 1s negligible; therefore’ the I,,r’/I,, and [,,p’/l,, terms 
are disregarded. 


In state space form, X’=AX+BU: 


pg) |%e Ye %_, sesO) pg ft % 
eu) U, Te, 
ete i, o Paar. 1, | 
"| IN, N, WN, 0 || in, Nu, Eo 
uM oa! tan( ) oo REF Lt O04] 


3. Longitudinal and Lateral-Directional EOM in State Space Form 


The full 8x8 state space equation is shown below: X’=AX+BU. The longitudinal and 
lateral-directional EOM are decoupled in this analysis. The longitudinal EOM need not be solved 


to analyze this problem. The longitudinal EOM are included for completeness. 
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4. Determination of the Non-Dimensional Stability Derivatives 


The non-dimensional force and moment coefficients were determined in wind tunnel 
testing and updated based on actual flight data. The equations and coefficients are found in the 
Operational Aerodynamic Data Book Volumes | and III [Ref. 12]. 


a. Longitudinal Stability Derivatives 


The Taylor series first order expansion of aerodynamic coefficients is given by: 
Cp: = Cro + Cpe + Cos(&)! 
Ch = Gren GrodGroloan 
Cr =O eee naa he aon 

The elevons are at trail (0°) at touchdown [Ref. 2]. Postlanding, the elevon is 
deflecting as required to balance the longitudinal moment equation and maintain a constant pitch 
attitude. However, elevon deflection is minimal until the start of derotation. At this time, the 
elevons ramp quickly downward as theta decreases at a rate of 1-2 deg/sec. Therefore, (6.)) is 
negligible in the two point stance and set to 0 in these three equations. 

The coefficients [Ref. 12] are defined as: 
Cor = Copasic FACphieh ar FAC real gas tAC (viscous intercon TOCOBE |r clevatorailern ia: anne 
FAC rpiteh rate FAC gig FAC GE FAC ysitts pod FAC (wehicle 
where the subscript ( ) denotes lift force (L), drag force (D), or pitching moment (m) and: 
Cypasic = Basic, full-scale, freestream, rigid body force or moment coefficient. 
AC high at = Change in force or moment coefficient due to high alt. effects (h>285,000 ft). 
AC real gas = Change in force or moment coefficient due to real gas effects. 
AC (yiscous interacuon = Change in force or moment coefficient due to viscous interaction (Vo. effects). 
ACopr = Change in force or moment coefficient due to body flap deflection. 
AC (levatorailcron = Change in force or moment coefficient due to elevator and aileron deflection. 
AC plunging = Change in force or moment coefficient due to rate of change of alpha. 
AC (pitch rate = Change in force or moment coefficient due to pitch rate. 
ACyiGc= Change in force or moment coefficient due to fully extended landing gear. 
AC vce = Change in force or moment coefficient due to proximity of the ground. 


AC silts pod = Change in force or moment coefficient due to addition of silt pods. 


AC vehicle = Change in force or moment coefficient due to geometry differences between vehicles. 

AC onigh att ANd AC (reat gas are not applicable below an altitude of 265,000 ft. and M 
< 6.0, respectively. AC viscous interaction 18 not a factor since V.. is zero. AC (ctevatoriaileron 1S aSSuMed to 
be zero since elevator/aileron deflection is considered negligible. AC piunging aNd AC (pitch rate are 
not a factor with the Orbiter on the ground in the two point stance since the pitch rate and angle of 
attack are essential constant until the initiation of derotation. This analysis will assume Orbiter 
Vehicle (OV) 103 which does not have a silts pod; therefore’ AC (sins poa 18 not applicable. The 
above equation can be reduced to: 
Cor = Copasict AC opr tAC gig tAC ce t+ AC yehicte 

Each of the terms in the above equation was obtained by linear interpolation of 
data extracted from a table at a = 10°, M=0.4 and a = 10°, M=0.25 to obtain data applicable for 
M=0.325 (353 fps) for both the lift and drag coeffiecients. Cza was determined by breaking each 
of Copasic. ACopr, ACouc, ACyce, and AC vehicle terms into Czo and Cy, components and summing, 
the C,,, contnibutions. Cp, was determined in the same manner. 


Since U is considered constant in this analysis, Cp and and C, are zero. 
Cy and C, are zero since the pitch rate is zero. The lift and drag coefficients pertinent to this 
g g 
analysis are summarize in Table 3. 


Table 3. Lift and Drag Coefficients 


Lift and Drag Coefficients Value (per °) 


0.54055 
0.05948 


0.00992 





b. Lateral-Directional Stability Derivatives 


The Taylor series first order expansion of the lateral-directional aerodynamic 
coefficients are provided for completeness. The Taylor series first order expansion of 


aerodynamic coefficients is given by: 


Cy = Cre =f Cyppi co Cylon a Cys(&)i 
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Cn = Clo + CipBi + Ctia{ da) + Cra( 6): 
Cm = Cno + Cupp + Créal dz)1 + Cus(&)1 

The lateral-directional bias coefficients are assumed to be zero for a symmetrical 
airframe: Cyo = Clo = Cno= 0. At touchdown, a software flag denoted ‘flat turn’ is set [Ref. 13]. 
When this flag is set, lateral motion is accomplished without elevon motion, using rudder only. 
Therefore, aileron deflection is negligible in the two point stance and (dz): are set to zero in the 
above three equations. The coefficients [Ref. 12] are defined as: 
Cor = Cosidestip AC Haiteron FC (nidder TAC (roll rate POAC aw males © Oaletpedia Sea tenele 
where the subscript ( ) denotes side force (y), yawing moment (n), or rolling moment (1), and: 
C(sidestip = Basic, full-scale, freestream force or moment coefficient due to sideslip angle. 
AC (aileron = Change in force or moment coefficient due to aileron deflection. 
AC (rudder = Change in force or moment coefficient due to rudder deflection. 
AC (roll rate = Change in force or moment coefficient due to roll rate. 
AC yaw rate = Change in force or moment coefficient due to yaw rate. 
AC (sits pea = Change in force or moment coefficient due to addition of silt pods. 
AC wehicle = Change in force or moment coefficient due to geometry differences between vehicles. 
AC (aileron May be dropped from this equation since aileron deflection is negligible. 

Each of the lateral-directional stability and control derivatives was obtained by 
linear interpolation of data extracted from a table at a@ = 10°, M=0.4 and a = 10°, M=0.25 to 
obtain data applicable for M=0.325 (353 fps). The side force, rolling moment, and yawing 


moment coefficients pertinent to this analysis are summarize in Table 4. 
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Table 4. Side Force, Rolling Moment, and Yawing Moment Coefficients 















Side Force, Rolling Moment, 


Value (per °) 
Coefficients 


and Yawing Moment 
| 
. eg 
ee 
. 0.0032 
| 


-0.2925 
: -0.00054 


+ 






| 


ee 


& 





LY 


tf 


eee | ee ae 


3 





of 





V. PROPOSED FCS MODIFICATIONS TO IMPROVE 


DIRECTIONAL CONTROL IN THE TWO POINT STANCE 


The rudder FCS susceptibility to PIO results from the slow surface rate, slow response 
time from rudder pedal input to rudder deflection, the power of the rudder, and the placement of 
the rudder with respect to the gear. Increasing the surface rate requires hardware modification. 
Analysis of hardware modifications are beyond the scope of this thesis. Transport delay might be 
decreased by shortening time delays in the FCS, linearizing the rudder pedal command shaper, 
and/or disabling Ny feedback. The power of the rudder might be reduced by decreasing the 
rudder pedal gain. 

Options to improve directional control may be categorized into expensive and 
inexpensive changes. Inexpensive options include those changes which require only software I- 
load or K-load change or a new piloting technique. Software I-loads and K-loads are constants 
which are not coded in the software and can be easily changed with no actual software code 
change. Expensive options include hardware changes, software changes which require a code 
change, and changes which require FCS recertification. 

Table 5 summarizes the problems with the current FCS, potential changes to the FCS to 
reduce susceptibility to PIO, the relative cost of the change, and whether the change is evaluated in 
this thesis. 


Table 5. Potential Changes to the FCS 


Slow rudder rate Hardware modification ae 


Slow rudder response | Decrease time delays Inexpensive 

























Linearize rudder command shaper 


Disable Ny feedback Expensive 


Powerful rudder Decrease pedal gain Inexpensive Re | 
Position of rudder Redesign Orbiter | Veryexpensive |No | 


Inexpensive 
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r% LINEARIZE THE RUDDER PEDAL COMMAND SHAPER 


The rudder command shaper was described in Section III.A.1.b. Recall that the shaping 


schedule is parabolic and conforms to the following function: 
DRMS = (0.042) DRMAN* + (0.131)DRMAN . 

The shaper can be linearized by changing 0.042 to 0 and 0.131. to 1.0. The above 
equation becomes: DRMS = DRMAN . This change provides a one-to-one response up to the 


limits of + 22.5° which may result in a more predictable rudder response. 
B. REDUCE THE RUDDER PEDAL GAIN 


The rudder pedal gain was described in Section III.A.1.c. Reducing the gain will reduce 
the g-command authority on the rudder and may reduce over-controlling of the Orbiter. The 
current gain value is 0.08574 for a touchdown speed of 207 KEAS. The proposal is to limit the 
gain, GNYDRM to 0.04 at 683 fps TAS and below. This change will reduce the gain by 47% 


which may reduce the over-controlling tendency. 
C DECREASE TIME LAGS 


The time lag between rudder pedal input and rudder deflection can be reduced by 
changing the Ny feedback time constant, described in Section III.A.1.f, from t = 0.8 s to t= 0.2 s. 
This change should result in quicker rudder response from the FCS which should reduce the 
susceptibility to PIO. 


D. DISABLE LATERAL ACCELERATION FEEDBACK 


Another option is to disable Ny feedback in combination with changing the Ny feedback 
time constant from T=0.8 s to t=0.2 s. This change would result in a more direct rudder FCS by 
allowing the rudder command to pass directly through the system with minimal delay. This 


change would require software code change and is relatively expensive. 
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VI. YAW FCS SIMULATION 


Directional control was evaluated by incorporating Orbiter lateral-directional contro! laws 
and longitudinal and lateral-directional EOM in MATLAB 5.0 and Simulink. The script program, 
‘eom.m', defines the initial conditions, mass properties, stability and control derivatives, and 
longitudinal and lateral-directional EOM in state space form. The code for 'eom.m' is contained in 
Appendix A. The code for the Simulink simulation is contained in Appendix B. The plotting 
program, 'FCSplot.m’, 1s contained in Appendix C. 


A. SIMULATION INITIAL CONDITIONS 


The case investigated is comparable to the STS-62 touchdown flight conditions. The 
commander input approximately +15° (right) rudder pedal deflection followed by -15° (left) 
rudder pedal deflection in the two point stance in an attempt to capture runway centerline. 
Landing speed was approximately 207 KEAS with minimal crosswind (< 2 knots). Table 6 


contains the initial conditions used tn the simulation. 


Table 6. Initial Conditions 


a 
ee 
ig Sure 
a 
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B. FCS CONFIGURATIONS EVALUATED 


Simulations were run to compare six FCS configurations as defined in Table 7. 


Table 7. FCS Configurations Evaluated 


[0 | Caren confguation sie) 















Reduce rudder pedal gain (rudgain) 
Cases 3,4,5 combined (done) 


C. SIMULATION 


The simulations were run in the time domain for 30 seconds using Dormand-Prince 
integrator. The simulated rudder pedal deflection included a rudder doublet of +15° rudder at 2 


seconds followed by -15° as depicted in Figure 19. The master Simulink block diagram is shown 


in Figures 20. 


Rudder Pedal Deflection 
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Figure 19. Rudder Pedal Deflection 
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Figure 20. Simulink Block Diagram of Rudder Doublet Effect on Directional Control 


The lateral-directional FCS is a subsystem of the entire system shown in Figure 20. The 
FCS is labeled "YAW FCS' in Figure 21. This subsystem contains the simulation code and block 
diagrams for the FCS. This block diagram and its internal block diagrams are shown in Figure 21. 
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Figure 21. YAW FCS Block Diagram 
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Vil. RESULTS 


Traces of Spcdai, 5;, Ny, B, r, p, yw and versus time are shown in Figures 22-28. All six 
configurations are shown on each plot for comparison of the proposed changes. All Simulink 
traces shown in Figures 22-28 are in response to the 15° doublet. The Simulink model was 
validated by comparing the results with actual Shuttle flight data and Shuttle Entry Simulator 
(SES) simulation results. 


A. MODEL VALIDATION 


To establish validity of the Simulink model, actual STS-62 flight data and SES simulation 
results were compared to the Simulink results. The STS-62 flight data and SES results are 
contained in Appendix D. The STS-62 data is compared to the Simulink ‘Baseline’ traces, since 
the onboard FCS is in the same configuration. The STS-62 CDR input a comparable 15° rudder 
pedal doublet in the two point stance. An SES evaluation, using the same initial conditions as this 
simulation and the STS-62 flight data, was run for comparison. Two SES runs were made. The 
SES 'Old FCS' traces contain the same FCS configuration as the Simulink ‘Baseline’ run. The SES 
'CF FCS' (Combined Fix) traces contain the same FCS changes as the Simulink ‘Optimum fix’. 

Rudder deflection traces for all three data sets (Simulink, STS-62 flight data, and SES) 
show similar magnitude of nearly 15° deflection, similar response time, and deflection rate. Ny 
traces exhibit non-minimal phase response as indicated by the initial low magnitude acceleration 
in the negative direction. Lag times are similar and the magnitude of the Ny traces are 
comparable. A companson of the sideslip angle traces indicate the same initial order-of- 
magnitude B resulting from the doublet. However, the Simulink trace shows an apparent long 
term directional instability which may be the reason why there is no B response in the opposite 
direction. The long term directional instability is reported to exist in the actual Orbiter in the two 
point stance. This instability is most likely not as strong as indicated in the Simulink model. 

Simulink yaw rate and roll rate traces also indicate the same magnitude and lag time as 
SES results and actual flight data. The SES ‘Old FCS' roll rate data is 'noisy' due to the gear 
model. 

This comparison of flight data and SES results with these Simulink traces results in a high 
level of confidence that this model is sufficient to analyze and compare the proposed FCS 


configurations for short term Orbiter directional response. 
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B. RUDDER DEFLECTION 


The rudder deflection plot in Figure 22 indicates that decreasing the time constant results 
in the quickest rudder response. Simply reducing Tt to 0, as in the ‘No Ny feedback’ option, yields 
a rudder response very similar to the pedal command. This option appears the best with respect to 
rudder surface control, yet the absence of Ny feedback generates other long term anomalies in 
yaw rate and siseslip stability. 

Reducing t to 0.2 with Ny feedback, as in the 'Shorter time constant’ trace, yields an 
acceptable rudder surface response with the benefits of Ny feedback. This trace shows 
approximately one half second of time delay by the end of the four second doublet. This time 
delay could also be considered as ~45° of phase lag. While 45° of phase lag is not desired, this 
option appears to be much better than some of the other options. The baseline configuration 
shows nearly two seconds of time delay before the rudder returns through the zero deflection 
position. This delay equates to nearly 180° of phase lag by the end of the maneuver. This 
response is obviously an unacceptable control response to a low frequency command. 

Linearizing the rudder pedal shaper, with no additional changes, would result in an 
increase in rudder deflection of 25% over the baseline. This larger rudder deflection for the 
telatively small rudder pedal command would aggravate the already difficult fine directional 
control. Linearizing the rudder with no other changes is not recommended. 

Reducing the rudder pedal gain, with no additional changes, is effective in decreasing the 
magnitude of the rudder surface response. However, the rise time in the deflection is too slow. 
The slow rise time distorts the predictability of the rudder response. This change does not effect 
the excessive phase lag. 

The ‘optimum fix’, which includes the t=0.2 time constant, linearized rudder shaper, and 
reduced rudder gain, result in the best shaped rudder response while retaining the Ny feedback. 


The deflection is shaped more like the rudder pedal input than all other Ny feedback cases. 
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Figure 22. Rudder Surface Deflection 


c. LATERAL ACCELERATION AT THE COMMANDER'S STATION 


The lateral acceleration plot shown in Figure 23 shows that eliminating Ny feedback 
results in greater Ny gain and more time required to return to a zero Ny condition. This option 
displays the worst long term directional stability characteristics. 

Both the 'Shorter time constant’ and ‘Linearized rudder’ fixes also result in an increase in g 


by 25% over the baseline FCS. This increase in lateral acceleration control power could lead to 
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directional oversensitivity. The ‘Shorter time constant" fix results in the least phase lag, while 
‘Linearized rudder’ yields the longest time delay of the Ny feedback options. ‘Reduced rudder 
gain’ results in an appropriate decrease in Ny gain. However, this change alone results in 


excessive time delay. 


The ‘Optimum fix’ has the least time lag and acceptable Ny gain characteristics. 


Lateral Acceleration at Commander Station 
025 





Q.2 


Ny (9) 


Figure 23. Lateral Acceleration at the Commander's Station 
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D. SIDESLIP 


The sideslip plot shown in Figure 24 shows the total rudder control power in the classic 
sense of degrees of sideslip per degree of rudder pedal. Most modem aircraft exhibit rudder 
control power on the order of 1:4, or one degree of sideslip for four degrees of rudder pedal. The 
‘Baseline’ Orbiter trace depicts rudder control power of 1:2.3. This apparent oversensitivity can 
lead to difficulty in precisely controlling sideslip. 

The 'No Ny feedback’ option exhibits excessive rudder control power with an additional 
long term instability. This long term response is unacceptable. The change that optimizes the 
classic 1:4 rudder control power ratio is the ‘Optimum fix’. 

The sideship plot also reveals the inherent long term directional instability of the Orbiter in 
the two point stance. This instability is due to the change in the center of rotation from the CG to 


the main wheels as discussed earlier. 
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Figure 24. Sideslip Angle 


E. YAW RATE AND HEADING 


Yaw rate and heading traces are shown in Figures 25 and 26. Heading was derived by 
integrating yaw rate. Long term directional instability is evident in both yaw rate and heading 
traces. The ‘Shorter time constant’ and ‘Optimum fix' options indicate the least phase lag of all 


options with Ny feedback. 
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Figure 25. Yaw Rate 
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Figure 26. Heading Angle 


ROLL RATE AND ANGLE OF BANK 

Adverse roll is evident as shown in Figures 27 and 28. This response is undesirable 
because adverse roll contributes to the apparent reduced directional control. The 'Reduced rudder 
gain’ results in the least adverse roll. However, this option results in significant phase lag. The 


“Optinum fix' offers reduce adverse roll and the least phase lag. The long term positive dihedral 


is evident in the long term positive roll rate. This roll rate 1s in response to the long term negative 
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Figure 27. Roll Rate 
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Figure 28. Angle of Bank 


G. NATURAL FREQUENCY AND DAMPING RATIO 


The so called ‘directional natural frequency’, Op, and damping ratio, Cp, were estimated 
from the plots after completion of the doublet. The damping ratio was estimated using the 


following equation: 


54 





Estimations from the plots indicate that @p = 0.6 rad/sec and Cp = 0.45. As expected, the 
natural frequency and damping ratio are low. This low natural frequency and low damping ratio 


can result in the tendency of the FCS to be susceptible to PIO. 


> 





VIII. CONCLUSIONS 


Orbiter landings indicate both long term directional instability and the potential for pilot 
induced oscillations during landing and rollout before nosewheel touchdown with the baseline 
(current) FCS. The Orbiter’s FCS requires improvements to mcrease directional contro! in the 
two point stance. This thesis described the control deficiency, potential improvements to the FCS, 
and evaluated six FCS configurations to determine their effect on directional control. The 
evaluation was performed by modeling the Orbiter’s postlanding lateral-directional contro! laws 
using MATLAB 5.0 and Simulink. The longitudinal equations of motion were include so that this 


imodel can more readily be adapted to evaluate longitudinal control problems. 
A. RECOMMENDATIONS 


Of the six FCS configurations evaluated, the option which provides the most 
improvement over the baseline configuration is the ‘Optimum fix’. This fix includes reducing the 
time lag t from 0.8 sec to 0.2 sec, linearizing the rudder pedal shaper, and reducing the rudder 
gain. This change minimizes the apparent phase lag, minimizes the adverse roll, reduces 
directional oversensitivity, and provides the most traditional rudder control power. This option is 
relatively inexpensive since only software K-load and I-load changes are required. Software code 
and hardware modifications are not required. The data indicates that PIO tendencies would be 
reduced over the baseline system if this option were implemented. 

This simulation was limited to analyzing the response of the FCS to a rudder doublet 
since the initial Orbiter directional difficulties were in response to a CDR doublet. The doublet 
input, in essence, is a frequency domain input. Therefore, frequency domain analysis would be 
appropriate. However, the validity of the Simulink model can more accurately be determined by 
comparing the time domain response traces of the individual states to actual flight data. Although 
this method was used as a diagnostic tool to study a frequency-concerned phenomenon, phase lag 
can be estimated from the traces. In addition, these changes can be comapared and evaluated with 
a pilot in the loop in time-based simulators such as the SES and Ames VMS. 

This author recommends further investigation of the ‘Optimum fix’ in the SES and Ames 
VMS. Frequency domain analysis would also be appropriate. In addition, the 
MATLAB/Simulink program used for this evaluation should be upgraded to include a gear model, 


Sy 


drag chute model, and wind model. The roll and pitch channels should be incorporated for 


completeness. 
B. USE OF COTS SOFTWARE FOR ENGINEERING EVALUATION 


This author concludes that commercially available engineering software packages such as 
MATLAB may be used to analysis proposed FCS changes prior to full up simulation in more 
expensive simulators. Future programs, such as X38, as well as Orbiter upgrade projects can 


benefit from the use of PC-based simulation. 


58 


APPENDIX A. MATLAB CODE FOR SCRIPT PROGRAM 


This appendix contains the Matlab code for the scnpt program, "eom.m". EOM.M 
defines the initial conditions, mass properties, stability and control derivatives, and longitudinal 


and lateral/directional EOM in state space form. 


$Initial conditions and EOM 


clear 
=— 2090; Ssurface area of the wing in ft%*2 
5-39.56; $chord length in ft 


e—-16.05/7; Swing span in feet 
W=228237; $weight in lbs 
g=32.174; $gravitational constant 
M=W/g; mass in slugs 
Txx=858677; 

77 -6869907.5; 

TzZz=7133448.5; 


sohift the moments of inertia from the cg to center of rotation 
(between the main gear) 
Pemg=iz2z74+M* ((1172.6-1084.15) /12)%2 


@b/h=0.05 at alpha = 10 


Mu=0.03; érunway rolling Mu 
Wile 2353; eer leveloclty an the x direction ft/sec 


rho= 0.002377; %sea level pressure 
theta=9*pi/180 %Pitch angle in radians 
alpha=10*pi/180 %Angle of attack in radians 
g=0.5*rho*U1*2 %Dynamic pressure at sea level 


sNon-Dimensional Stability and Control Derivatives 
CLa=0.05948; 

CDa=0.009924; 

CL1=0.54055; 

GDiI—OL16712; 

CLu=0; 

CDu=0; 


@yb—-0 .01705; 
Cyp=0; 

Cyr=0; 
Cydeltar=0.00138; 
enis——O 200205; 
Cip——-025189; 
Cilr=0.1890; 
Glideltar—0 00032; 
Cnb=0.0010; 
Cnp=0.028; 
Cnr=-0.2925; 
Cndeltar=-0.00054; 


eLongitudinal Dimensional Stability and Control Derivatives 


Xu=~g*S* (CDu+2*CD1) / (M*U1) 
Xa=-q*S* (CDat+2*CL1) /M 


a7 


Xxq=0; 
xXdeltae=0; 


24U=-G*S* (CLU+2* Chive 
Za=-q*S* (Chat2-Chija 
Zq=0; 

Zdeltae=0; 

Mu=0; 

Ma=0; 

Mg=0; 

Mdeltae=0; 


sLateral-directional Dimensional Stability and Control Derivatives 
Yb=q* S*Cyb/M 


Yo=-@" SB eyo (25 Uae S$1sS zero 
Yroq* > 8) Cys) Sis zero 
Ydeltaa=0; 


Ydeltar=q*S*Cydeltar/M 


Lb=g* S*B*Clb/Ixx 
Ep=q*S* BA 2* Clip (2 
LY=q*S7*B2* Clee ica 
Ldeltaa=0; 
Ldeltar=q*S*B*Cldeltar/Ixx 


Nb=q*S*~B*=Cnb/ lzmg 
Ne-O7S78 2 Cony (2 iznGg 
Nr=g*S*8°2*enr/ (2 -a7n0 aly 
Ndeltaa=0; 
Ndeltar=q*S*B*Cndeltar/Izmg 


Definition of dx/dt= Ax+Bu 


oxdot=[udot alphadot qdot thetadot betadot pdot rdot phidot] 
$x=[u alpha q theta beta p r phi] 
SU=[{DELTAE DELTAA DELTAR] 


A = [Xu Xa Xq -g*cos(theta) 0 0 
Zu/U1 Za/U1l Zq/U1+1 -g*sin(theta 
MusMalMa7 0" 0) Os0m 0; 

0. 0F00220'- 

O- Yb/Uls Yp/Ul Yr/U1—=1 6 seostahicia) 7 UE 
0 hb Epes. 0; 

O Nips Nie 10); 

00 Tetant{ theta) 76) 


O 0; 
1/U1 "0 6 Ose 


Cee @ 
2 Oya 
SS Set 


B=[(Xdeltae O 0; 
Zdeltae/Ul O 0; 

Mdeltae O 0; 

GO: -0; 

Ydeltaa/Ul Ydeltar/U1; 
Ldeitaa Ldeltar; 
Ndeltaa Ndeltar; 

O 0] 


Ce (a 


8 dy/dt = Cx+Du; C= identity matrix D=0 so that the output 
sof the state space is X 
Pi levyene) | 
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D=[ zeros (8, 3) ] 


B=[Xdeltae O 0; 
Zdeltae/Ul O 0; 
Mdeltae O 0; 
O70 0; 


Weerdeltaa/Ul Ydeltar/Ul; 


O Ldeltaa Ldeltar; 
O Ndeltaa Ndeltar; 
Omer 0] 
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APPENDIX B. MATLAB CODE FOR SIMULATION 


The code for the actual simulation, "baseline.m" is contained in this Appendix. The code 


is for the baseline systein. 


Model { 
Name “baseline” 
Version 2.00 
SimParamPage Solver 
SampleTimeColors off 
WideVectorLines off 
PaperOnentation landscape 
StartTime [Or 
StopTime US0s 
Solver ode45 
RelTol Wile 
AbsTol "le-6" 
Refine oe 
MaxStep 05) 
InitialStep "auto" 
FixedStep "auto" 
MaxOrder > 
OutputOption RefineOutputTimes 
OutputTimes "ty" 
LoadExtermalinput off 
ExternallInput "[t, uj" 
SaveTime off 
TimeSaveName “tout” 
SaveState off 
StateSaveName “xout" 
SaveOutput off 
OutputSaveName "yout" 
LoadinitialState off 
InitialState "xInitial" 


SaveFinalState off 
FinalStateName "xFinal" 


LimitMaxRows off 
MaxRows "1000" 
Decimation oe 
AlgebraicLoopMsg warning 
MinStepSizeMsg warning 


UnconnectedInputMsg = waming 
UnconnectedOutputMsg warning 


UnconnectedLineMsg warning 
ConsistencyChecking off 
ZeroCross on 
BlockDefaults { 
Onentation nght 
ForegroundColor black 
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BackgroundColor 
DropShadow 
NamePlacement 
FontName 

FontSize 

FontWeight 

FontAngle 

ShowName 

} 

AnnotationDefaults { 
HorizontalAlignment 
VerticalAlignment 
ForegroundColor 
BackgroundColor 
DropShadow 
FontName 
FontSize 
FontWeight 
FontAngle 

} 

System { 

Name 

Location 

Open on 

ScreenColor 

Block { 

BlockType 

Name 

Position 

Onentation left 

} 

Block { 
BlockType 
Name 
Position 
DeleteFcn 
PostSaveFcn 
Location 

j 

Block { 
BlockType 
Name 
Position 
Outputs 
q 

s 

Block { 
BlockType 
Name 
Position 
Onentation left 
Inputs 


white 

off 

normal 
“Helvetica” 
10 

normal 
normal 

on 


center 
middle 
black 
white 

off 
"Helvetica" 
10 

normal 
normal 


“baseline” 


(0, 3, 644, 447] 
white 


Derivative 
"BETADOT" 
{385, 320, 415, 340] 


Clock 

"Clock" 

[30, 330, 50, 350] 

“simclock BlockIsBeingDestroyed" 


"simclock Save" 


[12, 31, 119, 76] 


Demux 

“"Demux" 

[370, 154, 405, 236] 
"S° 


Mux 
"Mux" 
(310, 253, 335, 347] 


va 


} 

Block { 
BlockType 
Name 
Position 
Inputs 

} 

Block { 
BlockType 
Name 
Position 
Onientation 
NamePlacement 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
Save ToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Onientation 
Expr 

} 

Block { 
BlockType 
Name 
Position 
ExternalReset none 
InitialConditionSource 
InitialCondition 
LimitOutput off 
UpperSaturationLimit 
LowerSaturationLimit 
Show SaturationPort 
ShowStatePort 
AbsoluteTolerance 


down 


" 1 " 
off 
a @ hee 


left 


Mux 

"Mux 1" 

1235, 178. 270, 212] 
eye 


Scope 
"Ny" 
[215, 395, 245, 425] 


alternate 

off 

(188, 235, 512, 474] 
off 


on 
“auto” 

v2S " 

nS" 

off 
"ScopeData" 
on 

" 5000" 


Fen 
"Ny_" 
(240, 277, 275, 323] 


"(353 *u[3]}+353*u[2]+uf 1]*60)/2" 


Integrator 
"ps" 
[505, 380, 525, 400] 


internal 
" 0" 


"inf" 


" -inf" 


off 


off 
" auto " 
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Block { 
BlockType 
Name 
Position 
ShowPortLabels 
MaskType 
MaskDescription 
MaskHelp 


MaskPromptString 


MaskStyleStnng 
Mask Vaniables 


MaskInitialization 


MaskDisplay 


MaskIconFrame 


MaskIconOpaque 


MaskIconRotate 
MaskIconUnits 


Mask ValueStnng 


System { 
Name 
Location 
Open 


ScreenColor 


Block { 


BlockType 


Name 
Position 
Value 


j 
Block { 


BlockType 


Name 
Position 
Amplitude 
Frequency 
Phase 


SampleTime 


} 
Block { 


BlockT ype 


Name 


SubSystem 

"Pulse\nGenerator 1" 

[235435, 55,05] 

on 

“Pulse Generator" 

"Pulse Generator" 

"Pulse generator which ensures pulse transitions " 
"are hit. Provides a vector of pulses when " 

"the height is entered as a vector. Unmask " 

"to see how it works." 

"Period (secs):|Duty cycle (% of period):|Amplitu" 
"de:|Start time:" 

"edit,edit,edit edit" 

"period=@ | ;duty=@2;amplitude=@3;start=@4;" 


Nine 
"plot(0,0, 100, 100,[95,90,90,75,75,60,60,45,45,30," 


"30,15,15,5],[25,25,75,75,25,25,75,75,25,25,75,75" 
" 25,25]); \n " 


" \n " 


on 

on 

off 

autoscale 
"1000}.2]15|5" 


"Pulse\nGenerator1" 
[205, 368, 544, 559] 

off 

white 


Constant 
"Constant 1" 

[40, 135, 60, 155] 
"Oe 


Sin 

"Sine Wave" 

(35, 83, 65, 107] 

" ] " 

"2*p1./period" 

"-2* pi*(start./period + duty/200 - 1/4)" 


We ] " 


Step 
" Step " 
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Position 
Time 
Before 
After 


Block { 


} 


BlockType 


Name 
Position 


Threshold 


Block { 


} 


BlockType 


Name 
Position 
Port 


ieee. 05, 53] 
"start" 


“Oe 


“amplitude” 


Switch 

"Switch" 

PF 20565, 155,125} 
"sin(pi*(0.5 - duty/100))" 


Outport 

"out! " 
[250, 88, 280, 102] 
"4 l " 


Output WhenDisabled held 


InitialOutput 


Line { 


} 


SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


Line { 


j 


SrcBlock 
SrcPort 
DstBlock 
DstPort 


Line { 


j 


SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


Line { 


j 
} 


SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


Block { 
BlockType 


Name 


"Oo 


"Sine Wave" 
l 
[5, 0; 20, 0] 
"Switch" 
2 


"Switch" 


"out 1" 


"Constant 1” 
] 

[25, 0; 0, -30] 
"Switch" 
3 


"Step" 
(20, 0; 0, 35] 
"Switch" 


SubSystem 
"Pulse\nGenerator2” 
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Position 
ShowPortLabels 
MaskType 
MaskDescription 
MaskHelp 


MaskPromptStning 


MaskStyleStnng 
Mask Vaniables 
MaskInitialization ne 


MaskDisplay "plot(0,0, 


[25, 95, 55, 125] 

on 

“Pulse Generator" 

"Pulse Generator" 

"Pulse generator which ensures pulse transitions " 
“are hit. Provides a vector of pulses when " 

"the height is entered as a vector. Unmask " 

“to see how it works." 

"Period (secs):|Duty cycle (% of period):|Amplitu" 
“de:|Start time:" 

“edit, edit,edit,edit" 

"period=@ 1 ;duty=@2;amplitude=@3;start=@4;" 


100, 100,[95,90,90,75,75,60,60,45,45,30," 
"30, 15,1535),[25,25, 15, 75,25,29. 1 ee 
",25,25]); \n : 


"t \n " 
MaskIconFrame on 
MaskIconOpaque on 
MaskIconRotate off 
MaskIconUnits autoscale 
Mask ValueString "1000).2|-15|7" 
System { 
Name "Pulse\nGenerator2" 
Location [205, 368, 544, 559] 
Open off 
ScreenColor white 
Block { 
BlockType Constant 
Name "Constant 1" 
Position [40, 135, 60, 155] 
Value "0" 
} 
Block { 
BlockType Sin 
Name "Sine Wave" 
Position [35, 83, 65, 107] 
Amplitude ae 
Frequency "2* p1./period" 
Phase "-2*pi*(start./period + duty/200 - 1/4)" 
SampleTime "1" 
} 
Block { 
BlockType Step 
Name "Step" 
Position [53427;,65.535] 
Time "start" 
Before moe 
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After “amplitude” 
} 
Block { 
BlockType Switch 
Name "Switch" 
Position (1205 65, 155, 125] 
Threshold "sin(ni*(0.5 - duty/100))" 
} 
Block { 
BlockType Outport 
Name "out I" 
Position [250, 88, 280, 102] 
Port a 
Output WhenDisabled held 
InitialOutput SQ 
} 
Line § 
SrcBlock "Step" 
SrcPort | 
Points [20, 0; 0, 35] 
DstBlock "Switch" 
DstPort l 
} 
Line { 
SrcBlock "Constant 1" 
SrcPort 1 
Points [25, 0; 0, -30] 
DstBlock "Switch" 
DstPort 3 
} 
Line { 
SrcBlock "Switch" 
SrcPort ] 
DstBlock "out 1" 
DstPort ] 
} 
Line { 
SrcBlock "Sine Wave" 
SrcPort ] 
Points [5, 0; 20, 0] 
DstBlock "Switch" 
DstPort 2, 
} 
e 
; 
Block { 
BlockType Derivative 
Name "RDOT" 
Position [385, 260, 415, 280] 
Orientation left 
} 
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Block { 


BlockType Scope 
Name "Scope" 
Position (395, 435, 425, 465] 
Floating off 
Location [115, 431, 439, 670] 
Open off 
Grid on 
TickLabels on 
ZoomMode on 
TimeRange "auto" 
YMin "-5" 
Y Max ie 
SaveToWorkspace off 
SaveName "ScopeData" 
LimitMaxRows on 
MaxRows "5000" 
Decimation an 
SamplelInput off 
SampleTime "0" 
j 
Block { 
BlockType Scope 
Name "Scopel" 
Position (390, 380, 420, 410] 
Floating off 
Location [100, 161, 424, 400] 
Open off 
Grid on 
TickLabels on 
ZoomMode on 
TimeRange "auto" 
YMin "-5" 
YMax shy 
SaveToWorkspace off 
SaveName "ScopeData" 
LimitMaxRows on 
MaxRows "5000" 
Decimation he 
SamplelInput off 
SampleTime AOh 
j 
Block { 
BlockType StateSpace 
Name "State-space" 
Position [295, 180, 355, 210] 
A "A" 
B "B" 
C eC. 
D "Dp" 
XO "[0;0;0;0;0;0;0;0]" 
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} 

Block { 
BlockType 
Name 
Position 
Inputs 

} 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
VanableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Orientation 
VariableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
VaniableName 
Buffer 
Decimation 


oie 


" ] " 


ae 


n=l t 


left 


le 


re 1 9 


s1" 


aoe 1 " 


at (ta 


Sum 

" Sum" 
{100, 77, 130, 108] 
Wyott 


ToWorkspace 

“To Workspace" 
(95, 145, 155, 175] 
"doublet" 

"inf" 


ToWorkspace 

"To Workspace!" 
[175, 250, 235, 280] 
Sar 

"inf" 


ToWorkspace 
“To Workspace2" 
19585056155, 395] 


"t Ny " 
"inf" 


ToWorkspace 

"To Workspace3" 
[540, 110, 600, 140] 
“beta” 

int. 


ToWorkspace 
"To Workspace4" 
1570; 150,030,180] 


Pp 
"inf" 
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SampleTime 

j 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

j 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

j 

Block { 
BlockType 
Name 
Position 
VariableName 
Buffer 
Decimation 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 


ShowPortLabels 


System { 
Name 
Location 
Open 


ScreenColor 


Block { 


BlockType 


Name 


Bet 


" l " 
ve] ih] 


a lags 


oe 1 " 


ee 


al " 


oe 


eel " 


ToWorkspace 
"To Workspace5" 
[585, 220, 645, 250] 


I 
"inf" 


ToWorkspace 

"To Workspace6" 
[585, 335, 645, 365] 
” psi " 

"inf" 


ToWorkspace 

“To Workspace7" 
[70, 320, 130, 350] 
"time" 


"inf" 


ToWorkspace 

"To Workspace8" 
[585, 280, 645, 310] 
"aob" 

"inf" 


SubSystem 

"YAW FCS" 

[175, 179, 205, 231} 
off 


"YAW ECS" 


(4, 44, 648, 440] 


off 
white 


Inport 
in| ] " 


Position 
Port 

Port Width 
SampleTime 

j 

Block { 
BlockType 
Name 
Position 
Port 
PortWidth 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Port 
PortWidth 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Onientation 
ShowPortLabels off 
System { 


{35, 180, 50, 200] 
" l Ut 
vel “" 


ae ] a 


Inport 

Mine 2" 

[100, 290, 120, 310] 
Oe 

vei 


we ] 9 


[290, 270, 310, 290] 
a 
ue ] " 


=| " 


SubSystem 
"“COSALF" 
(320, 70, 350, 120] 
down 


Name "COSALF" 
Location [318, 77, 607, 247] 
Open off 
ScreenColor white 
Block { 
BlockType MATLABFcn 
Name "Cos(Alpha)" 
Position P5795, 133, 125} 
Orientation down 
MATLABFcn "cos" 
Output Width i 
} 
Block { 
BlockType Constant 
Name "alpha" 
Position [S8, 50, 132, 65] 
Onentation down 
Value "10*pi/180" 
} 
Block { 
BlockType Outport 
Name "out_1" 
Position (85, 150, 105, 170] 
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Onientation down 


Port role 
Output WhenDisabled 
InitialOutput iO 
} 
Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 
j 
Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


j 
} 

} 

Block { 
BlockType 
Name 
Position 
Value 

j 

Block { 
BlockType 
Name 
Position 
ShowPortLabels off 
System { 


held 


"Cos(Alpha)" 
] 


“out_ |} " 
| 


"alpha" 

] 
"Cos(Alpha)" 
J 


Constant 
"Constant" 

(60, 125, 80, 145] 
" ; 04" 


SubSystem 
"DRERR" 
(335, 255, 3652305] 


Name "DRERR" 
Location (180, 294, 476, 476] 
Open off 
ScreenColor white 
Block { 
BlockType Inport 
Name ‘ine 
Position (15, 75, 35, 95] 
Port We 
PortWidth "4" 
SampleTime "1" 
j 
Block { 
BlockType DiscreteFilter 
Name "YAW_RATE BBF" 
Position [60, 56, 245, 114] 
Numerator "[.7718479 -.3633895 .4980655]" 


Denominator "{1 -.3633895 .2699135]" 


SampleTime ae 
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Block { 
BlockType 
Name 
Position 
Port le 
OutputWhenDisabled 
ImitiaJOutput Oa 


Outport 
“out | " 
[265, 75, 285, 95] 


held 


} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

} 


} 
Block { 


BlockType 
Name 
Position 
Value 

} 

Block { 
BlockType 
Name 
Position 
Inputs 

} 

Block { 
BlockType 
Name 
Position 
Inputs 


} 

Block { 
BlockType 
Name 
Position 
Onrentation 
Value 

j 

Block { 
BlockType 
Name 
Position 


es "f 


oe 


"YAW RATE BBF" 


| 
Pout! " 
l 


1] in |] " 
l 


"YAW _RATE BBF" 


l 


Constant 

"DRT" 

[135, 245, 150, 265] 
100 


Product 
"GAIN" 
[135, 175, 160, 195] 


Product 
"GAIN_" 
(355, 185, 380, 205] 


Constant 

"GDRC" 

[475, 300, 515, 330] 
up 

"6.325" 


Constant 
"GERYFWDS" 
[485, 105, 505, 125] 


Wes; 


Onientation down 
Value poe 
j 
Block { 
BlockType SubSystem 
Name "GNYDRM_COMP" 
Position [105, 65, 135, 115] 
Onentation down 
ShowPortLabels off 
System { 
Name "“GNYDRM_COMP" 
Location [355, 132, 620, 292] 
Open off 
ScreenColor white 
Block { 
BlockType Lookup 
Name "GNY DRM Comp" 
Position (75, 43, 150, 107] 
Orientation left 
Input Values "[0 250 900 10000)" 
Output Values "{0.1 0.1 0.01 0.01)" 
} 
Block { 
BlockType Constant 
Name TAS: 
Position [170, 75, 190, 95] 
Orientation left 
Value "353" 
} 
Block { 
BlockType Outport 
Name "“GNYDRM" 
Position [15, 70, 35, 90} 
Orientation left 
Port sb 
OutputWhenDisabled held 
InitialOutput aor 
} 
Line { 
SrcBlock "“GNY DRM Comp" 
SrcPort ] 
DstBlock "GNYDRM" 
DstPort ] 
} 
Line { 
SrcBlock AS! 
SrcPort ] 
DstBlock "“GNY DRM Comp" 
DstPort I 


j 
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j 
Block { 


BlockType SubSystem 
Name “NY COMP" 
Position (2958170, 325, 220] 
ShowPortLabels off 
System { 
Name “NYCOMP" 
Location [227, 109, 594, 284] 
Open off 
ScreenColor white 
Block { 
BlockType Inport 
Name "in I" 
Position [15, 80, 35, 100] 
Port oh 
PortWidth vo1" 
SampleTime ".]" 
j 
Block { 
BlockType Gain 
Name "GRAY" 
Position (215, 77, 240, 103] 
Gain “sy 
j 
Block { 
BlockType TransferFcn 
Name “Transfer Fcn" 
Position BiSOmG5. 1/5. 1 ho! 
Numerator [hh 
Denominator "(0.8 1]" 
} 
Block { 
BlockType Outport 
Name “out_ 1" 
Position (260, 80, 280, 100] 
Port relies 
OutputWhenDisabled _held 
InitialOutput 0 
} 
Line { 
SrcBlock eine | 
SrcPort ] 
DstBlock “Transfer Fon" 
DstPort l 
} 
Line { 
SrcBlock "GRAY" 
SrcPort l 
DstBlock “out_1" 
DstPort ] 
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} 
Line { 


SrcBlock "Transfer Fcn" 
SrcPort l 
DstBlock "GRAY" 
DstPort l 
} 
j 
j 
Block { 
BlockType Product 
Name "Product 1" 
Position (STO. 1775 539213) 
Inputs 8 
j 
Block { 
BlockType Saturate 
Name "RUDCMD" 
Position [555, 190, 580, 210] 
UpperLimit elie 
LowerLimit "-27.1" 
} 
Block { 
BlockType SubSystem 
Name "RUDSHAPE" 
Position (80, 165, 110, 215] 
ShowPortLabels off 
System { 
Name "RUDSHAPE" 
Location [75, 132, 474, 272] 
Open off 
ScreenColor white 
Block { 
BlockType Inport 
Name cinele 
Position [15, 70, 35, 90] 
Port a 
PortWidth "1" 
SampleTime "-]" 
Block { 
BlockType Saturate 
Name "DR_ Limiter" 
Position (210, 66, 235, 94} 
FontSize 8 
UpperLimit pie 
LowerLimit veZe eS” 
} 
Block { 
BlockType Fen 
Name "DRman(0.42|DRman|+0. 131)" 
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Position 
Expr 
t 


5 
Block { 
BlockType 
Name 
Position 
Port i 
Output WhenDisabled 
InitialOutput uO? 
} 
Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 
} 
Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 
} 
Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


j 
j 


j 
Block { 


BlockType 
Name 

Position 
Floating 
Location 

Open 

Gnd 
TickLabels 
ZoomMode 
TimeRange 
YMin 

Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 


[105, 70, 145, 90] 


"(abs(u)*0.042+.131)*u" 


Outport 
" out_ 1 w 
(265, 70, 285, 90] 


held 


"DR_ Limiter" 
] 

“out_ 1" 

l 


"in 1 " 

1 

"DRman(0.42|DRman|+0. 131)" 
1 


"DRman(0.42/DRman|+0. 131)" 
l 

"DR_ Limiter" 

1 


Scope 

"Scope" 

[400, 130, 430, 160] 
off 

(188, 365, 512, 604] 
off 

on 

on 

yonly 

"auto" 

mgm 

oi 

off 

"ScopeData" 

on 

"5000" 

" ] " 

off 

"on 


i 


j 

Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Gnid 
TickLabels 
ZoomMode 
TimeRange 
Y Min 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid 
TickLabels 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid 


Scope 

"Scopel" 
[275, 110, 305, 140] 
off 
(188, 365, 512, 604] 
off 

on 

on 

on 

"auto" 

ng 

iS i] 

off 

"ScopeData" 

on 

"5000" 

Li | ” 

off 

non 


Scope 

"Scope2" 

[425, 300, 455, 330] 
off 

(188, 365, 512, 604] 
off 

on 

on 

on 

"auto" 

MES " 

"5m 

off 

"ScopeData" 

on 

"5000" 

nye 

off 

"oO" 


Scope 

"Scope3” 

[215, 115, 245, 145] 
off 

[188, 365, 512, 604] 
off 


on 
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TickLabels 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Gnid 
TickLabels 
ZoomMode 
TimeRange 
YMin 
YMax 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation 
SampleInput 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Inputs 

} 

Block { 
BlockType 
Name 
Position 
Inputs na" 

} 

Block { 
BlockType 
Name 
Position 
Inputs 


ip 


JE 


on 
on 

“auto” 

nS 

"gn 

off 
"ScopeData" 
on 

"5000" 

1] ] " 

off 

"on 


Scope 

"Scope4" 

[240, 250, 270, 280] 
off 

[188, 365, 512, 604] 
off 

on 

on 

on 

"auto" 

cS ¢ 

ase 

off 

"ScopeData" 

on 

"5000" 

"ye 

off 

KO" 


Sum 
" Sum" 
[420, 190, 440, 210] 


Sum 
"Sum1" 
35,0775 255,208) 


Sum 
"Sum2" 
[180, 175, 200, 195] 
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j 


Block { 
BlockType SubSystem 
Name "YBB FILTER" 
Position [445, 15, 475, 65] 
ShowPortLabels off 
System { 
Name "YBB_ FILTER" 
Location [55, 69, 652, 308} 
Open off 
ScreenColor white 
Block { 
BlockType Inport 
Name "in 1" 
Position [155 20357110) 
Port ole 
PortWidth "1" 
SampleTime re 
j 
Block { 
BlockType DiscreteFilter 
Name "Filter 1" 
Position [60, 76, 265, 124] 
Numerator "[.699558 -.82766 +.5599439 }" 
Denominator "LY e082 73 45540710" 
SampleTime Pies 
j 
Block { 
BlockType DiscreteFilter 
Name "Filter2" 
Position [S258 79 575.1710 
Numerator "[.940726 -1.04423 801291 }" 
Denominator "[1 -.96303 .66162]" 
SampleTime a 
j 
Block { 
BlockType Outport 
Name "out_ 1” 
Position [550, 90, 570, 110] 
Port lke 
OutputWhenDisabled held 
InitialOutput 20G 
j 
Line { 
SrcBlock "Filter2" 
SrcPort l 
DstBlock South 
DstPort l 
} 
Line { 
SrcBlock iat 
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SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


j 
j 


} 
Block { 


BlockType 


Name 
Position 
Port 


Output WhenDisabled 


InitialOutput 


} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 


[50, 0] 


(95, 0] 


[10, 0} 


| 
"Filter 1" 
] 


"Filter1" 
] 
"Filter2" 
] 


Outport 

“out. | " 

[605, 190, 625, 210] 
" 1 " 

held 

nos 


" Sum" 
] 


"Product 1" 
2 


"in 2 
] 


"Sum1" 
2 


" in | " 

| 
"RUDSHAPE" 
l 


"DRT" 
] 


" Sum2" 


"RUDSHAPE" 
] 
" GAIN" 
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DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 


tN 


"GAIN" 


" Sum2" 


"NYCOMP" 


15, 0;0, 5) 
"GAIN_" 


"COSALF" 
[0, 65] 

"GAIN " 

1 


"Producti" 

] 
"RUDCMD" 
] 


"GERYFWDS" 
1 

"Product!" 

] 


"GDRC' 
"Product" 
3 


“"RUDCMD" 


i out _ l "4 


" GAIN _" 
] 
[0, 0] 
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DstBlock 
DstPort 

} 

Branch { 
DstBlock 
DstPort 


} 

} 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

DstBlock 
DstPort 
} 
Branch { 
DstBlock 
DstPort 


} 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

Points 
DstBlock 
DstPort 

} 

Branch { 
Points 
DstBlock 
DstPort 
} 

} 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

DstBlock 
DstPort 


j 
Branch { 


SS0mM 


"Scope" 
] 
"Sum!" 
] 
[0, 0] 
"NYCOMP" 
J 
"Scopel" 
J 
"in 3 " 
] 
"DRERR" 
] 
"DRERR" 
l 
(20, 0] 
[15, 0] 
" Sum" 
2 
[0, 35] 
"Scope2" 
] 
"Sum2" 
] 
[0, 0] 
"Sum1" 


| 
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Points 
DstBlock 
DstPort 


} 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

} 


t 


s 
Block { 
BlockType 
Name 
Position 
Onentation up 
NamePlacement 
Floating 
Location 
Open 
Grid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation le 
SampleInput off 
SampleTime "0" 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
Save ToWorkspace 
SaveName 


(0, -55] 


"Scope3" 


"GNYDRM_COMP" 
i 
"GAIN" 


Scope 
"alfa" 
[415, 25, 445, 55] 


alternate 

off 

[188, 365, 512, 604] 
off 


yonly 

"auto" 

"O. 1 5 Ty 

“0: 25 " 

off 
"ScopeData" 
on 

"5000" 


Scope 

"aob" 
[550, 300, 580, 330] 
off 

[188, 235, 512, 474] 
off 


on 

"auto" 

WS0): ] " 

O02) 

off 
"ScopeData" 
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LimitMaxRows 
MaxRows 
Decimation le 
Samplelnput off 
SampleTime UO? 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation ole 
SampleInput off 
SampleTime OF 
j 
Block { 
BlockType 
Name 
Position 
Value 
} 
Block { 
BlockType 
Name 
Position 
Value 
} 
Block { 
BlockType 
Name 
Position 
Orientation up 
NamePlacement 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 


on 
"5000" 


Scope 

"beta" 

[S505, 155, 535, 185] 
off 

[188, 235, 512, 474] 
off 


on 
"auto" 

"0.4" 

"Oe 

off 
"ScopeData" 
on 

"5000" 


Constant 

"deltaa" 

[190, 135, 210, 155] 
eur 


Constant 

"“deltae" 
[190, 95, 210, 115] 
"oO" 


Scope 
"doublet" 
[145, 30, 175, 60] 


alternate 
off 


[188, 365, 512, 604] 
off 


on 
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TimeRange 
Y Min 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation i 
SampleInput off 
SampleTime "OE 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
YMax 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation ae 
SamplelInput off 
SampleTime rs 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation pole 
SampleInput off 


" auto" 

vey " 

a e 

off 
"ScopeData" 
on 

"5000" 


Scope 

"dr" 
[240, 230, 270, 260] 
off 

[188, 365, 512, 604] 
off 


yonly 

"auto" 

vet " 

ae 

off 
"ScopeData" 
on 

"5000" 


Scope 

"psi " 
[565, 375, 595, 405] 
off 

[188, 235, 512, 474] 
off 


on 
"auto" 

"=a. 05 " 

“Oe 5 = 

off 
"ScopeData" 
on 

"S000" 
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SampleTime Ne 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
YMax 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation ae 
SampleInput off 
SampleTime "0" 
j 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
YMax 
Save ToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation a 
SamplelInput off 
SampleTime cor 
} 
Block { 
BlockType 
Name 
Position 
Floating 
Location 


Open 


Scope 


q 
(455, 55, 485, 85] 
off 
[188, 365, 512, 604} 
off 


yonly 

"auto" 

ues " 

us n 

off 
"ScopeData" 
on 

"5000" 


Scope 

“rollrate" 
[545, 195, 575, 225] 
off 

[188, 235, 609, 474] 
off 


on 

"auto" 
"-0.05" 
"0.01" 

off 
"ScopeData" 
on 

"5000" 


Scope 

"theta" 

[475, 100, 505, 130} 
off 

[188, 365, 512, 604] 
off 
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Gnid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation ei: 
SampleInput off 
SampleTime "oO" 
} 
Block { 
BlockType 
Name 
Position 
Orientation left 
Floating 
Location 
Open 
Gnd on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
YMax 
Save ToWorkspace 
SaveName 
LimitMaxRows 
MaxRows 
Decimation le 
Sampleinput off 
SampleTime aCe 
j 
Block { 
BlockType 
Name 
Position 
Floating 
Location 
Open 
Grid on 
TickLabels on 
ZoomMode 
TimeRange 
YMin 
Y Max 
SaveToWorkspace 
SaveName 


yonly 

"auto" 

| "9 

" l : 5 9 

off 
"ScopeData" 
on 

"5000" 


Scope 


u 
[375, 65, 405, 95] 


off 
[188, 365, 512, 604] 
off 


on 
"auto" 

{ts 1 00 " 

"OM 

off 
"ScopeData" 
on 

"5000" 


Scope 

"yawrate" 
1555; 245, 585, 275] 
off 
[188, 235, $12, 474] 
off 


on 

"auto" 

al 0 

"9 05" 

off 
"ScopeData" 
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LimitMaxRows 
MaxRows 
Decimation el 
SamplelInput off 
SampleTime aie 
} 
Line { 

SrcBlock 

SrcPort 

Points 

Branch { 
Points 
DstBlock 
DstPort 

} 

Branch { 
DstBlock 
DstPort 

j 

} 
Line { 

SrcBlock 

SrcPort 

Points 

Branch { 
Points 
DstBlock 
DstPort 

} 

Branch { 

Points 

Branch { 
Points 
DstBlock 
DstPort 

} 

Branch { 
DstBlock 
DstPort 

} 

} 

} 
Line { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 
Line { 

SrcBlock 


on 
"5000" 


"Demux" 
8 
[95, 0; 0, 85; 30, 0] 


[0, -20] 
"To Workspace8" 
I 


"aob" 


"Demux" 
5 
[50, 0] 


[35, 0; 0, 130] 


"BETADOT" 
1 


[0, -30; 25, 0] 


[0, -25: 40, 0] 


"To Workspace3" 
l 


"beta" 


"Demux" 

4 

[40, 0; 0, -75] 
"theta" 

l 


"Demux" 


9] 


SrcPort 

Points 

DstBlock 

DstPort 

; 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

DstBlock 

DstPort 

} 

Branch { 
Points 
Branch { 

DstBlock 
DstPort 

} 

Branch { 
Points 
DstBlock 
DstPort 


} 
j 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


3 
[0, -5; 30, 0] 


q 
] 
"Demux" 


[20, 0] 
"alfa" 


"Demux" 


[10, 0] 


"YAW FCS" 
l 
[5, 0] 


"Mux |" 
3 
(0, 40] 


"ar 
] 


[-50, 0; 0, 20] 


"To Workspace" 
I 


"Mux" 
] 
"Ny." 
I 


"Mux l "f 

I 
"State-space" 
I 


OZ 


Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


Line { 
SrcBlock 
SrcPort 
Points 
Dst Block 
DstPort 


Line { 

SrcBlock 

SrcPort 

Points 

Branch { 

Points 

Branch { 
Points 
DstBlock 
DstPort 

} 

Branch { 
Points 
Branch { 

DstBlock 
DstPort 

} 

Branch { 
Points 
DstBlock 
DstPort 

} 

} 

} 

Branch { 
Points 
DstBlock 
DstPort 

j 

Branch { 
Points 
Dst Block 
DstPort 

} 

Branch { 
Points 


"deltae" 
] 

0) 
"Mux! tI] 
] 


"deltaa” 


[0, 50] 
"Mux!" 


"Demux" 
i 
(2550) 


[55, 0] 
[-5, 0; 0, 170} 
"PSI" 
I 
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“yawrate” 


] 


[30, 0] 
"To Workspace5" 
1 


[15, 0; 0, 155; -290, 0] 
"YAW FCS" 
3 


(0, 50] 


"RDOT" 


[40, 0; 0, 80] 
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DstBlock 
DstPort 
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} 

Line { 
SrcBlock 
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DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
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DstBlock 
DstPort 

} 

Branch { 
DstBlock 
DstPort 

} 

} 

Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

DstBlock 
DstPort 

} 

Branch { 
Points 
DstBlock 
DstPort 

} 

j 

Line { 
SrcBlock 
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Line { 
SrcBlock 
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Line { 
SrcBlock 
SrcPort 
Points 
Branch { 

Points 
DstBlock 
DstPort 
} 
Branch { 
Points 
Branch { 
DstBlock 
DstPort 

j 

Branch { 
DstBlock 
DstPort 
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Line { 
SrcBlock 
SrcPort 
Points 
Branch { 
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DstBlock 
DstPort 
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Branch { 
Points 
Branch { 
Points 
DstBlock 
DstPort 

j 

Branch { 
Points 
DstBlock 
DstPort 
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Line { 

SrcBlock 
SrcPort 
Points 
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Nas 
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[-110, 0; 0, -95] 
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"To Workspace2" 
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"doublet" 
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j 

Branch { 
Points 
DstBlock 
DstPort 


} 
} 
Line { 

SrcBlock 
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Branch { 
DstBlock 
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Branch { 
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DstBlock 
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Line { 
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APPENDIX C. MATLAB CODE FOR PLOT PROGRAM 


plot(time,doublet,'b’) 
title("(Rudder Pedal Deflection’) 
xlabel(‘Time (sec)') 
ylabel(‘dpedal (deg)’) 

figure 


plot(time,dr,'-',time,dr1 ,'--',time, dr2,':', time, dr3,'x' time,dr4,'.', time,dr5,'-.') 

grid 

title("Rudder Deflection’) 

xlabel('Time (sec)') 

ylabel(‘dr (deg)’) 

legend(‘Baseline’,"No Ny feedback','Shorter time constant','Linearized rudder’,"Reduced rudder 
gain’,‘Optimum fix’, 1) 

figure 


plot(time, Ny,'-',time, Ny ],'--',time, Ny2,':',time, Ny3,'x', time, Ny4,".',time, NyS,'-.’) 

grid 

title(‘Lateral Acceleration at Commander Station’) 

xlabel('Time (sec)') 

ylabel(‘Ny (g)’) 

legend(‘Baseline’,"No Ny feedback','Shorter time constant','Linearized rudder’,"Reduced rudder 
gain’,'Optimum fix’, 1) 

figure 


plot(time, beta* 180/pi,'-',time, beta 1 * 1 80/pi,'-- 

' time, beta2* 1 80/pi,':',time,beta3 * 180/pi,'x' time, beta4* 1 80/pi,'.', time, betaS* 1 80/pi,'-.') 

grid 

title('Sideslip') 

xlabel('Time (sec)’) 

ylabel(‘Beta (deg)’) 

legend(‘Baseline’,,No Ny feedback','Shorter time constant',"Linearized rudder','/Reduced rudder 
gain’,'Optimum fix',4) 

figure 


plot(time,r* 180/pi,'-',time,r 1 * 180/pi,'-- 

‘ time,r2* 180/pi,':',time,r3 * 180/pi,'x',time,r4* 180/pi,'.',time,r5* 180/pi,'-.' 

end 

title("‘Yaw Rate’) 

xlabel('Time (sec)') 

ylabel('r (deg/sec)') 

legend(‘Baseline',.No Ny feedback','Shorter time constant',Linearized rudder’,'Reduced rudder 
gain’,’'Optimum fix’, 1) 

figure 


plot(time, p* 180/pi,'-',time, p 1 * 180/pi,'-- 

' time, p2* 1 80/pi,':',time, p3 * 1 80/pi,'x' time, p4* 1 80/pi,".’,time, p5* 1 80/pi,'-.) 
grid 

title(‘Roll Rate’) 


oy 


xlabel('Time (sec)’) 

ylabel('p (deg/sec)’) 

legend(‘Baseline’,"(No Ny feedback’,'Shorter time constant’,'Linearized rudder',,Reduced rudder 
gain’,’'Optimum fix’, 1) 

figure 


plot(time,psi* 1 80/pi,'-',time, psi I * 180/pi,'-- 

' time, psi2* 180/pi,':’,time, psi3 * 180/pi,'x’ time, psi4* 1 80/pi,’.’, time, psi5 * 180/pi,'-.’) 

grid 

title(‘Heading’) 

xlabel('Time (sec)’) 

ylabel(‘psi (deg)’) 

legend(‘Baseline’,"No Ny feedback’,'Shorter time constant','Linearized rudder'’,,Reduced rudder 
gain’,'Optimum fix’, 1) 

figure 


plot(time,aob* 1 80/pi,'-',time,aob 1 * 1 80/pi,'-- 

' time,aob2* 180/pi,":’,time,aob3* 180/pi,'x',time,aob4* 180/pi,'.’,time,aob5*180/pi,'-.' 

grid 

title’ Angle of Bank’) 

xlabel('Time (sec)’) 

ylabel(‘phi (deg)’) 

legend(‘Baseline’,"No Ny feedback','Shorter time constant’,"Linearized rudder',"Reduced rudder 
gain’,'Optimum fix',4) 
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APPENDIX D. STS-62 FLIGHT DATA AND SES SIMULATION 


RESULTS 


The first four plots contain the STS-62 data. The remaining six plots are the SES results. 
For all cases, Main Gear Touch Down (MGTD) occurs at 80.3 seconds with derotation at 91 
seconds. The X scale reflects time from the Approach and Land (A/L) transition, which normally 


occurs around 10,000 feet. 
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